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We investigate fractional Brownian motion with a microscopic random-matrix model and introduce 
a fractional Langevin equation. We use the latter to study both sub- and superdiffusion of a free 
particle coupled to a fractal heat bath. We further compare fractional Brownian motion with the 
fractal time process. The respective mean-square displacements of these two forms of anomalous 
diffusion exhibit the same power-law behavior. Here we show that their lowest moments are actually 
all identical, except the second moment of the velocity. This provides a simple criterion which enables 
to distinguish these two non-Markovian processes. 
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Diffusion is one of the basic non-equilibrium phenom- 
ena. Normal diffusion is well described in the theory of 
Brownian motion as a Gaussian process that is both local 
in space and in time. It is characterized by a mean- 
square displacement which is asymptotically linear in 
time, (a;^) = 2Dt, where D is the diffusion constant [Q. 
However, a growing number of experimental observations 
show that more complex processes, in which the mean- 
square displacement is not proportional to t, also occur in 
nature. Anomalous diffusion has for instance been seen 
in micelle systems [H, in two-dimensional rotating flows 
in porous glasses in actine networks but also 
on capillary surface waves Q, in strongly coupled dusty 
plasmas [Q, and more recently in intracellular transport 
Anomalous diffusion finds its dynamical origin in 
non-locality, either in space or in time. A well-known 
example of a process which is non-local in space is Levy 
stable motion, for which the mean-square displacement is 
actually infinite due to the occurrence of very long jumps 

. In this Letter we focus on processes which are non- 
local in time and whence show memory effects. Specifi- 
cally, we shall discuss and compare fractional Brownian 
motion (fBm) and the fractal time process (ftp) . 
These two forms of anomalous diffusion are fundamen- 
tally different (see below). Yet, they are difficult to tell 
apart experimentally, since both yield a mean-square dis- 
placement of the form (a;^) oc t", a ^ 1. It is for instance 
still an open question whether the long-range correla- 
tions observed in nucleotide sequences are to be 
interpreted in terms of fBm or ftp type DNA walks [p^ . 
In this paper we aim at providing a simple criterion which 
permits to distinguish between these two non-Markovian 
processes. 

The very difference between fBm and ftp is best il- 
lustrated by looking at their diffusion equations. The 
solution of the diffusion equation for fBm p6[ 

-P,^Ux,t) = aDt"-^—P,^^{x,t) (1) 

is easily found to be the Gaussian distribution 
PfB„(a;,t) = exp(-xV4i:)P)/[47ri:)P]i/2. FBm thus de- 



scribes Gaussian transport. It is important to note that 
Eq. (|l|) is thereby local in time (there is no memory ker- 
nel) . The non-Markovian character is expressed through 
a time-dependent diffusion constant, Da{t) = aDt"~^. 
In contradistinction, the diffusion equation for ftp |f 7| 
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contains a memory kernel and the distribution function 
Pitp{x,t) is hence non-Gaussian. The solution of (||) is 
given by Pftp[a;,z] = exp{-\x\z°'^^/D^/^)/[2D^/^z'^-"^'^] 
in Laplace z space. In time, P[t-p{x,t) has been ex- 
pressed in closed form in terms of a Fox function or 
a one sided Levy stable distribution By introduc- 

ing further the Riemann-Liouville fractional derivative 
(0<A<1) M 



f{r)dT 



(3) 



Equation (0) can be rewritten as a fractional diffusion 
equation pll 



Q_l_Petp{x, t) 
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Both Eqs. (|l|) and (^) reduce to the normal diffusion 
equation when a = l. 

We begin our discussion of fBm by introducing a frac- 
tional Langevin equation. It is worthwhile to point out 
that the Langevin and the phase-space descriptions of 
Brownian motion are no longer fully equivalent in the 
non-Markovian regime of interest here. As recently dis- 
cussed by Calzetta et al. , the Langevin equation con- 
tains more information and thus appears more fundamen- 
tal. We then apply this fractional Langevin equation to 
study in some detail the anomalous diffusion of a free 
particle coupled to a fractal heat bath. In particular, we 
evaluate the first two moments of both the position and 
the velocity of the particle, which we express in terms of 
Mittag-Leffler functions. Finally, we compare with the 
results obtained recently for ftp by Metzler and Klafter 
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for < a < 1 ||2g] and by Barkai and Silbey for 1 < a < 2 
psf by using a fractional Klein-Kramers equation. We 
find that fBm and ftp satisfy the same generalized Ein- 
stein relation. Moreover, their lowest moments are all 
equal, except the second moments of the velocity. 

We examine the dynamics of the Brownian particle 
with a microscopic random-matrix model. Random- 
matrix theory has already been successfully applied in 
the context of anomalous diffusion in Refs. p5|-|27[|. We 
thus consider a system S weakly coupled to a fractal heat 
bath B via a random-matrix interaction p7| , p8| . The 
coupling is chosen linear in the position x of the system. 
The generic form of the Hamiltonian is given by 



H = Hs(E>tB + '^s®Hb +x(E)V 



(5) 



where Hs — /2M + U{x) is the Hamiltonian of the sys- 
tem, Hb describes the bath and F is a centered Gaussian 
random band-matrix. It is assumed that initially the sys- 
tem and the bath are uncorrelated and that the latter is 
in thermal equilibrium at temperature (3 = {kT)~^ . The 
variance of the random interaction is further taken to 



have the form 27 



[p{ea)p{eb)\'- 
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Here e^'s denote the eigenenergies of the bath Hamilto- 
nian {HB\a) — £a\a)), Aq is the strength of the coupling, 
A the bandwidth and p{£) is the density of states of the 
bath, which is locally written as /9(e) = /3oexp(/3e). As 
shown in |^^, the variance (||) gives rise to subdiffusion 
when a < 1 and to superdiffusion when 1 < a < 2. The 
coupling to the bath is characterized by the bath correla- 
tion function which is defined as K{t) — {V{t)V{Q))B ~ 
K'{t) +iK"{t). Here V{t) = exp{iHBt)V exp{-iHBt) 
and {■ ■ ■)b denotes the thermal average. After perform- 
ing the average over the random-matrix ensemble, K{t) 
is found to be simply the Fourier transform of the vari- 
ance Vab^ with respect to £{,. In the following we con- 
sider the limit of high temperature and large bandwidth, 
1 ^ A <C kT. Using the variance ^ we then obtain 

W{t) = 2A„Tia)cosi^)t-\ K^{t)^^'^. (7) 

We see that the time dependence of K(t) follows an 
inverse power law. This presence of a long tail leads 
to long-time correlation effects in the dynamics of the 
Brownian system Note that for a — 1, the Fourier 

transform of (^) reads K'{t) = 2'nAQ5{t) and normal 
Brownian motion is recovered. Moreover, in the limit 
of weak coupling considered here, we have shown in H^] 
that the random band-matrix model can be mapped onto 
the oscillator bath model There is a one-to-one 

correspondence between the variance Vab^ of the ran- 
dom band-matrix model and the spectral density func- 
tion J{oj) of the oscillator bath model. The generalized 



Langevin equation corresponding to the random-matrix 
Hamiltonian (|^) can therefore be readily written as |^ 

Mi{t) + M [ 7(< - T)±{T)dT + U'{x) = ^{t) , (8) 
Jo 

where ^(t) is a Gaussian random force with mean zero 
and variance (C(i)C(O)) — K'{t), and j{t) is a damping 
kernel which obeys MkTj{t) = K'{t). This last relation 
is often referred to as the second fluctuation-dissipation 



theorem 31 1 . Remark that the Langevin equation is com- 
pletely determined by the real part K'{t) of the bath 
correlation function. Furthermore, in the limit of weak 
coupling, the dynamics described by the Hamiltonian (|^) 
is Gaussian and one can show that the corresponding dif- 
fusion equation is precisely given by Eq. (|l]). Using again 
the fractional derivative (j^), we may rewrite Eq. (g) in 
the form of a fractional Langevin equation. We obtain 
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Mi + Mia Q^^it) + = m , 



(9) 



where we have defined 7^ = 7rAo/3/(M sin(a7r/2)). The 
fractional Langevin equation describes both subdif- 
fusion for < a < 1 and superdiffusion for 1 < a < 2 [pll . 
As a simple application of the fractional equation M), 
we now concentrate on the free particle and accordingly 
set U{x) = 0. In this case, the solution of the Langevin 
equation is easily obtained by applying Laplace trans- 
form techniques |^ . We find 

x{t) =xo + voB^{t) + f B„it - r)C(T)dT , (10) 
Jo 

where {xo,vo) are the initial coordinates of the par- 
ticle and By{t) — J^Cv{t')dt' is the integral of the 
(normalized) velocity autocorrelation function C.u(t) = 
{v{t)v)/{v'^). The Laplace transform of Cv{t) is given by 



CJz] 



z + Jlz] 



Z + JaZ" 



(11) 



where ^\z] is the Laplace transform of the damping ker- 
nel. Eq. (|l l|) is known as the first fluctuation-dissipation 
theorem j^. By taking the inverse Laplace transform, 
the velocity autocorrelation function can be written as 



C^t) ^ E^-ai-lat^-") 



(12) 



Here we have introduced the Mittag-Leffler function 
Ea{t), which is defined by the series expansion |]3^ 



Ea{t) 



n— ^ ^ 



(13) 



The function Ea{t) reduces to the exponential when 
a — 1. The asymptotic behavior of the Mittag-Leffler 
function (O) for short and long times is respectively 
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given by ~ exp(t) and — (tr(l — a))~^ . For the ve- 
locity autocorrelation function (|lj) this yields a typical 
stretched exponential behavior at short times 



C'v (t) ~ exp ■ 



2-Q 



t«l/(7a)" 



r(3-a)' 

and an inverse power-law tail at long times 



^a^ia - 1) 



t > 
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The result ( p^ ) has already been derived in Ref. p7[ , 
where it has been shown to induce the "whip-back" ef- 
fect. After time integration, we finally get from Eq. (|l^ ) 

where we have used the generalized Mittag-Leffler func- 
tion Ea^f3(t) defined as Q 



r(an + P) 



(17) 



In the long-time limit, the generalized Mittag-Leffler 
function satisfies Ea^pit) ~ —{tr{f3 — a))^^. Accord- 
ingly, By (t) exhibits a decay of the form 



when < — > oo 



(18) 
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We emphasize that the solution (|l^) of the fractional 
Langevin equation in the force free case is completely 
specified by the knowledge of the function B^ (t) . 

Let us now turn to the evaluation of the lowest mo- 
ments of the position and the velocity of the free particle. 
The mean displacement and the mean-square displace- 
ment are readily deduced from Eq. (|lO|). We find 



(x) =Xo+VQtE2^a,2{~'lct'^ ") 



t^oo la r(a) 



(19) 



and 



2kT 



M 



2kT 



t->oo JaM r(l 



(20) 



In the last equation, thermal initial conditions have been 
assumed (xo — 0, Vq — kT/M). In addition, one may 
easily verify that (a;)^, (x^) and Cv{t) satisfy the general 
Green-Kubo relation 



(21) 



which is known from linear response theory pl| |. In a 
similar way, one can compute the first and second mo- 



ments of the velocity from the time derivative of Eq. ( 10 ) 
This results in 



and 



t^oo 7, 



^2-a(-7at'~") 

kT ( 



^0 ,„-2 



(22) 
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£;2-a(-7at'-")] } 
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We observe that (t;^) decays like (i"~^)^ for large t. A 
common remarkable property of the above calculated 
mean values is their slow relaxation towards equilibrium 
as given by the (generalized) Mittag-Leffier function. 
This has to be contrasted with normal Brownian mo- 
tion where all this quantities display an exponential de- 
cay. Let us now discuss the generalized Einstein relation 
which relates driven and free process . We consider a 
particle initially at rest (xg =Vq= 0) and seek the mean 
position {x)f as a function of an externally applied con- 
stant force U{x) — -~xF9{t). From the Langevin equa- 
tion we easily find 



d{x) 



dt 



(24) 



where the velocity autocorrelation function C„ (t) is given 
by Eq. ([l^). Equation (^) together with the Green- 
Kubo relation ( ^l|) for {x^)o in the force free case, then 
yields the generalized Einstein relation for fBm 



F 
2kT 



(25) 



It is interesting to note that the validity of the Einstein 
relation ( p5|) has been recently verified experimentally 

11,13. 

We now come to the comparison of fBm with ftp. 
Barkai and Silbey have investigated superdiffusive ftp 
with a fractional Klein-Kramers equation they inferred 
from a generalized Rayleigh model . For the free par- 
ticle, a direct comparison p8| between their results and 
our Eqs. (|l9|)-(p5|) shows that the mean displacement 
(|l9|), the mean-square displacement (20), the velocity's 
first moment ( ^^ and the velocity autocorrelation func- 
tion ([l^ ) are identical for the two processes. This means 
in particular that fBm and ftp satisfy the same Green- 
Kubo relation (pT[). Moreover, both fBm and ftp obey the 
same generalized Einstein relation (^5|). Although fBm 
and ftp are fundamentally different processes, we thus 
notice that they share strikingly common features. How- 
ever, the second moments of the velocity are different. 
For convenience, we quote their equation (2.18) which 
reads (in our notation) 



(1-2),.^ =Z;2^2-a(-27„i2-") 



kT 



|{l-^2-a(-27ai'~")} 



(26) 
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We see that for ftp, the second moment of the velocity re- 
laxes asymptotically like This is in sharp contrast 
to the fBm result Eq. ( p3| ) which exhibits a much faster 
decay. It is also worth noting that Eqs. ( p2| ) and ( |2^ ) 
reduce to the same (exponential) expression for a = 1. 
On the other hand, subdiffusive ftp has been studied by 
Metzler and Klafter by using a fractional Klein-Kramers 
equation derived from a non-Markovian generalization 
of the Chapman-Kolmogorov equation. A comparison 
with their results for the force free case leads to similar 
conclusions as in the superdiffusive regime. Many exper- 
iments on anomalous diffusion have measured either the 
mean-square displacement or the generalized Ein- 
stein relation [^0. However, the latter do not allow 
to distinguish fBm and ftp, as we have just shown. In 
contrast, the variance of the velocity offers a clear dis- 
tinction between these two processes as exemplified by 
Eqs. (23) and (p6|). This is the main result of this Letter. 

In summary, we have investigated fBm within a 
random-matrix approach and introduced a fractional 
Langevin equation which applies for both sub- and su- 
perdiffusion. We have studied the anomalous dynam- 
ics of a free particle coupled to a fractal heat bath 
and performed a comparison between fBm and ftp. We 
have found that these completely different forms of 
non-Markovian anomalous diffusion share many common 
characteristics. In particular, they satisfy the same gen- 
eralized Einstein relation and their lowest moments are 
all equal with the exception of the second moment of the 
velocity. 

We thank E. Barkai for reading the manuscript. 
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